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We present a systematic study of the shape function for inclusive _B-meson decays in 
the heavy-quark limit, which is based on the QCD equations of motion and heavy-quark 
symmetry, and takes into account the cusp divergence due to radiative corrections. 
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Inclusive _B-meson decays, such as semileptonic B Xuii" decays and penguin- 
induced B Xsj decays, are of special interest because they are sensitive probes 
of electroweak parameters as well as new physics. Especially for the precise determi- 
nation of the CKM matrix element \Vub\, the region near the kinematic endpoint for 
the lepton enerry spectrum, Eg ^ EJ"^^ = tub/'^.^ in the decay B Xu(-i> plays an 
important role™ avoiding large backgrounds from the decays into charmed particles, 
B Xcdv. In this endpoint region, the hadronic decay products evolving from the 
M-quark have large energy but small invariant mass, so that the hadronic final state 
becomes jet-like with coUinear interactions of an outgoing light quark. This implies 
that the corresponding differential decay rates are expressed as the light-cone expan- 
sion, analogously with deep inelastic lepton-nucleon scattering (DIS) cross sections: 
the leading term in the light-cone expansion is described by an analogue of the 
leading twist in the DIS, i.e., by a factorization formula where a structure function 
with all nonperturbative, long-distance ('^ 1/Aqcd) contribution is convoluted with 
the perturbatively calculable function. This fact orizat ion formula is valid in the 
leading power of AqcD/^^fc and to all orders in Q!s,EEEland the structure function, 
called the shape function, is expressed as the _B-meson matrix element of a bilocal 
light-cone operator in the heavy-quark effective theory (HQET). 

{B{v)\K{tn)¥e'<^^o<i^^-^^^''^KmB{v)) = f{t) = j doje'^'f{oj) , (1) 

where is the 4- velocity of the B-meson {v^ = 1), and denotes a light-like vector 
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{v? = 0, V ■ n ~ \), pointing in the direction of the final-state hght-quark jet. hy{x) 
denotes the effective 6-quark field, b{x) ~ exp(— imfjU • x)hv{x), satisfying '^h^ — 
f{t) and /(w) denote the shape functions in the coordinate and (residual) 
momentum representations, respectively. The shape function is process independent 
so that describes other inclusive B decays, e.g., the photon energy spectrum of 
B —t Xg^. For brevity, we do not show the path-ordered gauge factors connecting 
the constituent fields in the following. 

The bilocal operator of Q is renormalized at the scale /i. Recently it has been 
emphasized that the radiative corrections induce a Sudakov- ty pe str ong scale de- 
pendence on (^3), governed by the "cusp anomalous dimension" pl^l^l and, as a re- 
sult, f{t) is singular for short distances, f{t) ~ — {a sCp /t^)^o^ {it ^) as t — > 0; 
for /(w), this implies that "radiative tail" is generat ed for uj <C — /x, and all non- 
negative moments Jdujuj^f(uj) become UV divergent This contamination of 
the short-distance singularity, that cannot be removed by the usual renormaliza- 
tion, suggests^ that the bilocal operator of has to be further decomposed as 

K{tn)hy{0) = Y.^o{t,^l)0{^l) (2) 
o 

with the Wilson coefficient functions Ca{t, /i) that represent the "hard components" 
in integrating out the modes associated with the mass scales > 1/t; they are 
perturbatively calculable as Co = C'^'' -I- {as/Tr)C^^ + ■ ■ and separate out all 
short-distance singularities of Cq ^I, ^) ^ CFlog^(ii/i) as t ^ 0. On the other 
hand, the local operators 0(/i) represent the "soft components" and contain all 
nonperturbative effects. The OPE (0) is useful for t/x < 1, and an independent set 
of O and the corresponding coefficients Co can be determined by matching both 
sides of (O at the scale /z ~ 1/t. We perform the matching, order by order in as, and 
derive the sophisticated operator structure of that was previously unknown. 

at LO are independent of /i; therefore, Cq'^ does not have any noncanonical 
^-dependence singular for t ^ 0. Thus, in the tree-level matching of (O, hy{tn)hy{0) 
can be Taylor expanded about t = 0, and one can treat the corresponding local 
operators JiyD^^ ■■■D^.hy. However, this actually requires a compHcated task to 
disentangle JiyD^^ ■ ■ ■ D^.hy by a basis of "canonical operators" using the HQET 
equations of motion (EOM), iv ■ Dhy = 0, and has been feasible for the first few 
operators.!^ Here we solve the EOM constraints for the nonlocal operators, and 
Taylor expand the results at the final step!^ in the exact operator identity, 

d - - ^ - 

v^- — hy(x)hy{0) = hy(x)v -D hy(0)+i / duu hy{x)gGn,,(ux)v^x''hy(0) (3) 
dxt^ Jo 

with G^i, — {i/g)[Df^,Di,], the first term in the RHS vanishes using the EOM. 
Because x'^ is not restricted on the light cone, the corresponding constraint equation 
for Xf^ — > tUf^ is derived by combining with the light-cone expansion, 

hyix)hy{0) = [hyix)hyiO)]^^ + ^ / " TT^ni " M ( O) + (^(x^) , (4) 

^ JO ^ oXfj^oxf^ 
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where all local operators arising in the Taylor expansion of the first term are traceless 
and the following terms represent the corresponding "trace part" . We get 

t^/i„(in)/i„(0) + /i„M/Ji'(0)-/Ji.(0)/i^(0) = J(t) , (5) 



J{t)=t^l du{2iuhy{tn)gG^i,{utn)v^n'^hy{0) — — — - — hy{ux)hy{0) 
Jo \ uox^oxf^ 



(6) 



Here the second and third terms in the LHS of jSJ give the analogue of Nachtmann's 
correction in the DIS, due to deviation from the light-cone in 

Disentangling the last term of © with another operator identity from the EOM for 
heavy-quark and gluon fieldsj^the "source" term © becomes (i?_L^ = D^^ — v^v-D) 



2 



J{t) = -J dT{Thy{Tn)'v'^^hy{{)) +*'^ y duu^hy{Tn)g'^t''q{uTn)t''jiq{uTn)hy{0) 

-I- 2T^y duu J dsshy{Tn)gGf,^{uTn)n''gG'^f{sTn)nphy{0)j , (7) 

where the first term is nonlocal version of the kinetic energy operator. Eqs. lO, © 
show that the t-dependence, i.e., the longitudinal- momentum dependence, of the 
shape function is controlled by higher dimensional operators, that represent Fermi 
motion effects and 4-particle correlations with additional quarks and gluons inside 
the B meson, and are immediately integrated and Taylor expanded about t = 0: 

1 /■* 

hy{tn)hy{0) = hy{0)hy{0) + - dx J [t) 

t Jo 

= hyhy - —hyD\hy " — ky g'^ t"^ qt"^ fiqky + — — hyD\hy 
6 3d 120 \ 

-jhygh-qf^jf^qD'^^hy - hygGp^n'^gG^^Puphy^ +■■■ , (8) 

where we have simplified the tensor structure of the operators in the RHS using 
Lorentz invariance; for the j-th power P , the local operators of dimension j -I- 3 con- 
tribute. The expression for the general j-th power showJ^that an increasing number 
of operators participate for increasing j (i.e., ~ operators for j ^ 1). The re- 
sult © gives |2Jl with an operator basis {O} = {hyhy, hyD^hy, hyg^f^qf^i/iqhy^ ■ ■ ■} 
and the corresponding LO coefficients C^-* . 

Utilizing this operator basis, we proceed to the one-loop matching of (jSJ to ob- 
tain Cq\ For the bilinear-type operators O — hyhy and hyD^hy, the calculation 
is straightforward. We consider the two-point heavy-quark Green function with the 
insertion of Q, and compute the one-loop corrections for the nonlocal and local 
operators corresponding to LHS and RHS of (jJl, respectively: for hy{tn)hy{Q) this 
was calculated in the coordinate representation in RefP^ the number density o;^ 
erator hyhy and the kinetic energy operator hyD\hy obey non-renormalization.'^ 
We get C£iji,/i) = -CF[\og^{itfie^^) - \og{ittxe^^) + 5^724], cj'j^^^^ (t, ^i) = 
— i^C^^-*^ (t, ^)/6 in the MS scheme. Going over to the a;-representation via 1^, this 
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result coincides with the corresponding coefficients obtained in Ref|2i, and is consis- 
tent with those calculated in the momentum cutoff schemePl Matching for the other 
operators, O ~ hvg^t'^qt°"^qhy, ■ ■ •, can be treated similarly, but the complications 
arise because of multiparticle nature of the relevant operators. 

Substituting all above results into (O, we get the OPE for the shape function 
valid up to O (a^) and O {{AQcut)^cts) corrections: for low dimensional operators 
O = hvhy,hyD\hy, the coefficients Co are obtained at the NLO accuracy, and 
matrix elements of those operators directly give the familiar HQET parameterJHI 
at the scale -IGeV, as {B{v)\hyhy\B{v)) = 1, {B{y)\KD\hy\B{v)) = -Ai. For 
higher dimensional operators O = hyg^t°^qt°^t/iqhy, hyD'^hy, hyg^t'^qt°-jfj^qD'^hy, 
hygGfj.i^n'^gG'^Pnphv, ■ ■ ■, the coefficients Co are obtained at the LO, and ma- 
trix elements of those operators are related to the nonperturbative parame- 
ters at IGeV scale through evolution governed by the corresponding anoma- 
lous dimension. There are some simple estimates of those nonperturbative 
parameters.Elxhe anomalous dimension of KgH'^qt'^M. is known at one-loopP 
while the anomalous dimension matrix describing the renormalization mixing among 
{hyD\hy^ hvg^t°'qt'^Ji^iqD'^hy, hygC fj^^n" gC^^nphy} is not known. We note that 
matrix elements of the j+3 dimensional operators generate O ((AqcdO"' ) power cor- 
rections in the OPE (0) because canonical size of the corresponding nonperturbative 
parameters is ^ J^qcD' regions with rrib » these power corrections are 

enhanced compared with the usual subleading corrections, which are suppressed by 
powers of Aqcd / fnt- Physically, this corresponds to the situation in which the decay 
spectra in the endpoint region is smeared over a range A with rnj, ^ A > Aqcd- 
When t c± 1/Aqcd, corresponding to A ~ Aqcd, the higher-order power corrections 
in (0) are not suppressed, and the OPE has to be resummed. Under certain assump- 
tions, our differential equation ^ gives a "resummed" solution,^ whose behavior 
is consistent with a modeP guided by IR renormalon analysis. 

In conclusion, our results for ^ show that the shape function is a much more 
complicated object than the simple momentum distribution of the 6-quark inside the 
_B-meson, due to novel behavior of perturbative as well as nonperturbative nature. 
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